Superradiant light scattering from a moving Bose-Einstein condensate 
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We investigate the interaction of a moving BEC with a far detuned laser beam. Superradiant 
Rayleigh scattering arises from the spontaneous formation of a matter-wave grating due to the 
interference of two wavepackets with different momenta. The system is described by the CARL- 
BEC model which is a generalization of the Gross-Pitaevskii model to include the self-consistent 
evolution of the scattered field. The experiment gives evidence of a damping of the matter-wave 
grating which depends on the initial velocity of the condensate. We describe this damping in terms 
of a phase-diffusion decoherence process, in good agreement with the experimental results. 
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The experimental realization of Bose-Einstein conden- 
sates (BECs) with alkali trapped atoms has opened the 
possibility of investigating several fundamental aspects 
of quantum mechanics in macroscopic, i.e. many parti- 
cle systems lj. In superradiant Rayleigh scattering the 
coherent nature of the condensate leads to strong corre- 
lations between successive scattering events, as shown in 
the pioneering work of Ketterle and coworkers Q ■ This 
process was then the basis for the first demonstration 
of phase-coherent matter wave amplification Q . The ef- 
fect studied by Ketterle is an example of a spontaneous 
formation of a regular density grating in an atomic sys- 
tem, arising from a collective instability as in the Col- 
lective Atomic Recoil Laser (CARL) |J. In the ab- 
sence of thermal broadening (as it happens in a BEC), 
CARL appears as a promising source of macroscopically 
entangled or number-squeezed atom-atom and/or atom- 
photon systems 0, 0,0- However, in a real BEC sev- 
eral effects due, for instance, to spontaneous emission, 
inhomogeneous broadening and collisions, may seriously 
inhibit the CARL process and destroy the coherence in 
the matter wave field §• The control of decoherence in 
the photon-BEC interaction would be a significant step 
toward the achievement of macroscopic entanglement of 
coherent matter waves. 

In this paper we investigate both theoretically and ex- 
perimentally the influence of the initial velocity of the 
condensate on superradiant Rayleigh scattering. In the 
experiment we produce an elongated BEC of rubidium 
atoms and expose it to a single off-resonant laser pulse di- 
rected along the condensate symmetry axis. The laser is 
far detuned from any atomic resonance and the only scat- 
tering mechanism present is Rayleigh scattering 2] . In an 
elongated condensate a preferential direction for the scat- 
tered photons emerges, causing superradiant Rayleigh 
scattering. In this regime the atoms, initially scattered 
randomly, interfere with the atoms in the original mo- 
mentum state creating a matter-wave grating with the 
right periodicity to further scatter the laser photons in 
the same mode. Both the matter-wave grating and the 
scattered light are then coherently amplified. In our ge- 



ometry photons are back-scattered with k s sa — k, where 
k is the wave-vector of the laser photon, and the atoms 
move away from the original condensate with a relative 
momentum 2Tik in the direction of the laser beam. The 
efficiency of the process is limited by the decoherence 
between the original and the recoiled atomic wavepack- 
ets causing the damping of the matter- wave grating. We 
identify two different mechanisms for decoherence, one 
resulting from Doppler and mean field broadening of the 
matter wave field 0, and the other due to phase dif- 
fusion. The latter mechanism, dependent on the energy 
separation between the initial and final states of the sys- 
tem 0, [nl, can be controlled by initially setting the 
condensate into motion. In particular, we observe that 
phase diffusion decoherence vanishes when the initial con- 
densate momentum is such that after the interaction with 
the laser beam the scattered atomic wavepacket has the 
same kinetic energy of the original condensate in the lab- 
oratory frame. 

The experiment is performed with a cigar-shaped con- 
densate of 87 Rb produced in a Ioffe-Pritchard magnetic 
trap by means of RF-induced evaporative cooling. The 
axial and radial frequencies of the trap are uj z /2ir = 
8.70(7) Hz and uj r /2ir = 90.1(4) Hz respectively, with the 
z-axis oriented horizontally. We tune the atomic velocity 
after the end of the evaporative ramp by inducing a col- 
lective dipole motion of the condensate along the z-axis. 
The dipole oscillation is excited by non-adiabatically dis- 
placing the center of the magnetic trap. When the con- 
densate has reached the maximum velocity in the mag- 
netic potential, we suddenly switch off the trap and let 
the cloud expand with a horizontal velocity proportional 
to the displacement of the trap (see Fig^). We apply 
a square pulse of light along the z-axis, 2 ms after the 
release of the condensate, when the magnetic field of the 
trap is completely switched off, and the atomic cloud 
has still an elongated shape. After 2 ms of expansion 
the radial and axial sizes of our condensates are typi- 
cally 10 and 70 /im, respectively. The pulse length is 
controlled with an acousto-optic modulator. The light 
comes from a diode laser red-detuned 13 GHz away from 
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FIG. 1: Schematics of the experimental procedure. The con- 
densate is set in motion by a sudden displacement of the mag- 
netic trap center (A). When the condensate reaches the de- 
sired momentum po we switch off the magnetic trap and flash 
the atoms with a far off resonance laser pulse directed along 
the condensate symmetry axis (B). After an expansion time 
allowing a complete separation of the momentum components 
(28 ms) we take an absorption image of the atoms (C). 



the rubidium D2 line at A = 780 nm and has an intensity 
of 1.35 W/cm 2 corresponding to a Rayleigh scattering 
rate of roughly 5xl0 2 s _1 . The linearly polarized laser 
beam is collimated and aligned along the z-axis of the 
condensate. In this geometry the superradiant light is 
backscattered and the self-amplified matter-wave propa- 
gates in the same direction of the incident light. In order 
to minimize spurious reflections we have aligned the laser 
beam at a nonzero angle with respect to the normal to 
the vacuum cell windows. After an expansion of 28 ms, 
when the two momentum components are spatially sep- 
arated, we take an absorption image of the cloud along 
the horizontal radial direction. 

In Fig.^3 we show a typical absorption image in which 
the left peak is the condensate in its original momentum 
state po and the right peak is formed by atoms recoil- 
ing after the superradiant scattering at po + 2hk. The 
spherical halo centered between the two density peaks 
is due to non-enhanced spontaneous processes, i.e. ran- 
dom isotropic emission following the absorption of one 
laser photon. From a 2D-fit of the pictures assuming a 
Thomas-Fermi density distribution we extract the num- 



ber of atoms in both the original and the recoiled peaks. 

To understand the observed behavior we analyze in 
detail the self-consistent interaction between a coher- 
ent electromagnetic wave of amplitude a and a coher- 
ent matter wave field , 3>. The evolution of the system is 
described by the following 1-D CARL-BEC model, i.e. 
a Gross-Pitaevskii model generalized to include the self- 
consistent evolution of the scattered radiation amplitude: 
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In Eqs.Q and Q, a(t) = (e V/2hw s ) 1 / 2 E s (t) is the di- 
mensionless electric field amplitude of the scattered beam 
with frequency u s , g = (Q/2A)(cjd 2 /2heoV) 1 ' 2 is the 
coupling constant, f2 is the Rabi frequency of the laser 
beam with a frequency uj detuned from the atomic reso- 
nance frequency loq by A = uj — ujq, d — e-d is the electric 
dipole moment of the atom along the polarization direc- 
tion e of the laser, V is the volume of the condensate, 
N is the total number of atoms in the condensate and 
5 = lu — uj s . The matter wave field is normalized such 
that J dz\^\ 2 = 1. The nonlinear term in Eq.Q can be 
neglected since the experiment has been performed after 
expansion. The last term in right-hand side of Eq.QJ 
represents the self-consistent optical wave grating, whose 
amplitude depends on time according to Eq. (J2J . The first 
term in the right-hand side of Eq. J2J represents the self- 
consistent matter- wave grating. Eq.J2J) has been written 
in the "mean-field" limit, which models the propagation 
effects with a damping term where k « c/2L and L is 
the condensate length . 

If the condensate is much longer than the ra- 
diation wavelength and approximately homogeneous, 
then periodic boundary conditions can be assumed 
and the wavefunction can be written as ^S(z,t) — 
E„ c n (t)u n (z)e~ mSt , where u n {z) = (2/A) 1 / 2 exp(2inkz) 
are the momentum eigenstates with eigenvalues p z — 
n(2hk). Introducing the density matrix p m ^ n = c m c* n 
and u n = Aujjin 2 — Sn, where lor = hk 2 /2m is the recoil 
frequency, we obtain from Eqs.Q) and (J2J: 
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The last term added in Eq. © describes a phase-diffusion decoherence process, whose amplitude is characterized 
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by a constant t. This term, fundamental to describe 
our experimental results, arises from a (5-correlated gaus- 
sian noise on the eigenenergies of the system and causes 
the decay of the off-diagonal matrix elements, so that 
the density matrix becomes diagonal in the basis of 
the recoil momentum states. Eq.Q may be written 
as a master equation j 1 ll ] for the density operator p = 

E m , n p™,™l m >H : 

dp i , ~ , t r * 
i = -^[H,p]-- 2 [H ,[Ho,p]}, 



(5) 



where H = H + V, H = Ahuo R p 2 - hSp, V = 
ihg(a*e 2lkz — h.c.) and p — p z /2hk is the normalized mo- 
mentum operator with eigenstates \n) and eigenvalues n. 
The phase destroying term with the double commuta- 
tor of the Lindblad form in the right-hand side of Eq. JSJ) 
generates the damping term added in Eq.JSJ). It has ap- 
peared in many models of decoherence and induces diffu- 
sion in variables that do not commute with the Hamilto- 
nian, preserving the number of atoms in the condensate. 
In this term we have neglected the interaction V in the 
weak-coupling limit g 2 N/n <C u>r. 

In our experimental conditions the superradiant 
Rayleigh scattering involves only neighboring momentum 
states, i.e. transitions from the initial momentum state 
Pa = n(2hk) to the final momentum state (n + l)2hk. 
In this limit, our system is equivalent to a two- level sys- 
tem and Eqs.© and reduce to a set of equations for- 
mally equivalent to the well-known Maxwell-Bloch sys- 
tem j! E3: 
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where S = p„,„ + ie jA " 4 , A = ae lA " 4 , W = P n - P n+1 is 
the population fraction difference between the two states 
(where P n = p nj „ and P n + P n+1 = 1), A„ = u - u s - 
4ujR(2n + 1) and the decoherence rate 7„ is given by: 
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To the decoherence rate 7„ we have added an ex- 
tra term 70 taking into account other coherence de- 
cay mechanisms, as for instance Doppler and inhomo- 
geneous broadenings of the two-photon Bragg resonance 
0, El- We note that in Eq.J&jl, S represents half of 
the amplitude of the matter- wave grating. In fact, if 
^ w c n u n (z) + c n+ iu„ + i(z), the longitudinal density is 
|*| 2 w (2/A){l + 2Re[S* exp(2ikz - iA n t)]}, which de- 
scribes a matter wave grating with a periodicity of half 
the laser wavelength. The main result is that the second 
term of Eq. © , arising from a phase diffusion decoherence 
mechanism, depends on the frequency detuning between 
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FIG. 2: Left) Time evolution of population in the original 
condensate (empty circles) and in the recoiled wavepacket 
(filled circles) for different pulse durations. The solid line is a 
fit with the hyperbolic tangent I1UH predicted by the theoreti- 
cal model, the dotted line is just one minus the fit curve. The 
momentum of the original condensate is set to —Tik (top), 
(center) and +hk (bottom). Right) Plot of the atomic density 
profile after interaction with a 250/iS pulse for the three cases 
of original momentum as on the left. The laser detuning and 
intensity are 13 GHz and 1.35 W/cm 2 respectively. 



the incident and scattered radiation beams and on the 
initial momentum of the condensate, po — n(2hk). We 
observe that the velocity-dependent term of the decoher- 
ence rate is invariant under Galilean transformation. In 
fact, in a frame moving with respect to the laboratory 
frame, the shift of po compensates the Doppler shift of 
the frequency difference u> — uj s . 

Our experimental conditions match those for the su- 
pcrfluorescent regime 1 1 3| . in which the field loss rate 
k is much larger than the coupling rate g^/N . In this 
regime, for t 3> we can perform the adiabatic elim- 
ination A ~ gNS/(n — iA n ). The analytical solution for 
the fraction of atoms with initial momentum pa = n(2hk) 
is 

"2^ G 

x {l + tanh[(G-2 7 „)(t-t )/2]}, (10) 

where G = 2g 2 Nn/(n 2 +Al) is the superradiant gain and 
to is a delay time. In our experiment k 3> A n , so that 
G w 2g 2 N/K hence independent from the atomic veloc- 
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FIG. 3: Decoherence rate as a function of the initial momen- 
tum of the condensate. The solid line is a fit of the experimen- 
tal data with a parabola centered in po = — hk, as expected 
from the theoretical model. 

ity. Ea. (|10f> assumes the threshold condition G > 2^ n i.e. 
the gain must be larger than the decoherence rate. Fig[21 
reports the population fractions of the initial wavepacket, 
P n , and of the scattered wavepacket, P n +i, as functions 
of the laser pulse duration for three different initial mo- 
menta po. ^From the fits we extract the values of G and 
7„ for different po. The measured value of G = 19(3) 
ms _1 does not appreciably depend on pq. On the con- 
trary we observe a strong dependence of the decoher- 
ence rate 7„ on the initial momentum pq. In FigJSJwe 
plot the experimental points for the decoherence rate 
7„ as a function of the initial momentum of the atoms. 
The data show a parabolic behavior in good agreement 
with the prediction of Eq.JSJ if one assumes uj = uj s 
in the laboratory frame, with fit values 7o=4.2(2) ms' 1 
and t=2.4(2)-1CP 7 s. The theoretical calculation of the 
linewidth of the Bragg resonance |16| for our experimen- 
tal parameters predicts a value 70 ~ 3 ms _1 close to the 
value obtained from the fit in FigEJ Notice that the de- 
coherence rate is minimized for pq — —hk. Indeed, if the 



initial momentum is — hk, after scattering the atoms have 
the same kinetic energy and, with the above assumption 
for the scattered light frequency u> s , the phase-diffusion 
decoherence term in Eq.@ is zero. This identifies a sub- 
space which is decoherence free with respect to the phase 
destroying process [l7j . 

In our experimental apparatus it is very difficult to 
completely exclude the presence of pump light diffused 
by the vacuum cell windows. We have measured the dif- 
fused light present in our chamber in the direction oppo- 
site to the pump beam to be 1.0(3) x 10~ 5 of the pump 
beam intensity. This amount of light is of the same order 
of magnitude of the equivalent input noise of the super- 
radiant process [18j|. This justifies the assumption that 
u) s = u> used to fit the experimental data of Fig. [3] We 
remark that this back diffused light can not explain our 
results in terms of Bragg scattering since the width of 
the Bragg resonance is one order of magnitude smaller 
than the range of momenta explored in our experiment 
and shown in Fig. [3] Furthermore the hyperbolic tan- 
gent dependence of the atomic population in Fig. [21 can 
only be explained by the self consistent amplification of 
the matter wave grating and of the backscattered light 
as described in Eqs.QJ and 

In conclusion, we have studied the superradiant light 
scattering from a moving Bose-Einstein condensate. The 
efficiency of the overall process is fundamentally limited 
by the decoherence between the two atomic momentum 
states. With the assumption for the scattered light fre- 
quency ui s = uj in the laboratory frame, Eq.@ predicts 
the parabolic behavior of the decoherence rate as a func- 
tion of the initial momentum, in agreement with the ex- 
perimental results as shown in Fig[3] The fully quantized 
version of the CARL-BEC model offers the possibility of 
investigating the realisation of macroscopic atom-atom 
or atom-photon entanglement 0, . The control of de- 
coherence obtained in this work represents a significant 
step in this direction. 
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